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ABSTRACT 


- — ---> Dispersed  surface  waves  that  arrive  at  a  seismograph 
by  traveling  other  than  a  great-circle  path,  or  have  been 
delayed  by  reflection  along  the  path,  interfere  with  the 

• 

first  arriving  great  circle  waves.  The  problem  of  finding 
the  correct  group  velocity  dispersion  in  the  presence  of 
such  multipath  interference  can  be  simplified  by  the  use  of 
phase  equalization  filters  based  on  assumed  dispersion  curves. 
Group  velocity  dispersion  curves  calculated  using  zero  cross¬ 
ing,  multiple-filter,  or  moving-window  methods  can  be  per¬ 
turbed  slightly  with  each  perturbation  yielding  a  different 
phase  equalization  filter.  Cross-correlation  of  the  '•chirp'* 
implied  by  the  filter  with  the  dispersed  wave  train  produces 

« 

a  time  function  whose  real  spectrum  can  be  computed.  The 
process  is  repeated  with  different  trial  ''chirps'*  until  the 
integral  of  the  real  part  of  the  spectrum  is  maximized.  The 
associated  dispersion  curve  is  then  considered  to  be  a  *best 
estimate1*  of  the  true  curve.  Application  of  the  technique  to 
synthetic  seismograms  has  shown  a  resolution  significantly 
greater  than  existing  methods  when  multipath  effects  are 
present.  The  use  of  this  method  is  demonstrated  with  tele- 
seismic  Rayleigh  waves  recorded  at  LASA.  'V 
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INTRODUCTION 

If  Rayleigh  waves  propagated  only  in  a  laterally 
homogeneous,  layered  half-space  the  dispersion  curves 
obtained  by  the  calculation  of  periods  and  velocities 
from  the  peak  and  trough  or  zero  crossing  method  (Ewing 
and  Press,  1952;  McDonald  et  al.,  1974)  or  from  Fourier 
analysis  (Sato,  1955;  Block  and  Hales,  1968;  Landisman 
et  al,,  1969;  Dziewonski  et  al . ,  1969)  would  be  the  true 
dispersion  curves.  However,  dispersion  data  obtained  by 
these  techniques  are  almost  invariably  discontinuous; 
many  show  non-least  time  arrivals  (energy  refracted  or 
reflected  from  lateral  boundaries)  which  interfere  with 
the  wave  that  has  traveled  a  true  great  circle  path. 

The  effect  of  superposition  of  multipathed  signals  with 
the  primary  signal  is  to  produce  a  spectrum  with  "holes" 
and  "buildups"  at  various  frequencies.  In  addition, 
significant  errors  are  introduced  into  the  estimation  of 
the  dispersion  curve.  In  order  to  obtain  a  correct 
dispersion  curve  for  a  particular  travel  path  for  either 
Rayleigh  or  Love  waves,  the  effects  of  multipathing  must 


be  removed. 
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CHIRP  FILTERING 

Basic  Principles 

The  concept  of  phase  equilization  (or  chirp) 
filtering  can  be  utilized  in  the  estimation  of  dispersion 
curves  for  Rayleigh  waves  which  have  been  complicated  by 
multipathing .  If  an  initial  estimate  of  the  unknown  curve 
is  made  (e.g.  the  world-wide  average  for  continental  or 
oceanic  paths  of  Oliver  (1962)  would  suffice),  a 
corresponding  chirp  wave  form  can  be  generated  for  this 
estimate  by  the  method  of  McDonald  et  al .  (1974).  Dis¬ 
persion  curves  used  in  this  process  should  be  continuous, 
monotonically  increasing  functions?  thus  the  resulting 
chirps  will  be  smoothly  dispersed  signals  of  constant 
amplitude  with  phase  determined  by  the  estimate  of  the 
dispersion  curve. 

Let  c(t)  be  the  chirp  derived  from  an  estimated 
dispersion  curve  and  g(t)  the  recorded  signal.  Then 

G(t)^C(t)  »  jfg(t)  c(t+T)dt  =  H(T) 


r  - 
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where  H(T)  is  the  cross-correlation  function.  The 
Fourier  transform  of  H(T)  can  be  written 

G(«.u)  C(;u)  exp  i  [0g(>.o)  -  0c  (uj)J 
or 

G(to)  C  (tu)  |  cos  |0g(u/)  -  0c(<jj)J+  isin|0g(uj)  -  0c(u.')jj. 

Clearly  the  real  part  of  the  transform  of  H(T)  will  be 
maximized  when  the  phase  spectra  of  the  chirp  and  signal 
are  equal.  In  this  case  the  amplitude  spectrum  C(io)G(a>) 
will  be  real  and  even.  The  process  of  maximizing  the 
integral  of  the  re^l  spectrum  corresponds  in  the  time 
domain  to  maximizing  the  peak  value  of  the  correlation 
function.  Assuming  that  no  multipath  interference  is 
present,  this  maximization  can  be  achieved  in  practice 
by  varying  the  phase  of  the  chirp  until  it  matches  that 
of  the  signal.  The  "pseudoautocorrelation"  function 
(or  Paf)  is  then  defined  as  the  particular  H(T)  obtained 
when  phase  equalization  has  been  achieved.  A  Paf  has 
two  characteristics  which  set  it  apart  from  other  cross¬ 
correlation  functions:  (1)  it  has  the  maximum  peak 

possible  for  any  constant  amplitude  assigned  to  the 
chirp,  and  (2)  it  is  an  even  function  about  this  peak. 
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The  second  condition  results  because  the  frequency 
domain  equivalent  of  the  Paf  is  real  and  even,  and  the 
transform  of  a  real  and  even  function  is  also  real  and 
even  (Bracewell,  1965).  A  cross-correlation  function 
other  than  a  Paf  can  be  expressed  as  a  sum  of  its  even 
and  odd  parts — a  fact  which  accounts  for  its  peak  being 
smaller  in  magnitude. 

If  g(t)  contains  no  multipaths,  the  dispersion 
curve  could  be  determined  by  use  of  a  trial  and  error 
method  to  equalize  the  phase  spectra  0g(u-')  and  0c  (u>) . 

The  trial  dispersion  curve  whose  chirp  maximizes  the 
integral  of  the  real  part  of  the  transform  of  the  result¬ 
ing  cross-correlation  function  would  be  the  most  accurate 
estimate.  Of  course,  if  g(t)  has  no  multipaths  this 
analysis  is  unnecessary  since  the  dispersion  curve  could 
be  adequately  determined  by  other  means. 

In  the  case  where  g(t)  is  the  combination  of 
multiple  signals  arriving  at  different  times  it  can  be 
shown  that  H (T)  will  have  a  series  of  peaks  corresponding 
to  each  component  of  the  total  signal.  Consider,  for 
example,  the  case  where  g(t)  is  composed  of  a  true  great 
circle  path  signal  and  one  multipathed  arrival  delayed 
by  a  time  tQ. 
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Then 

g(t)  -  g^t)  +  Ag2(t-t^)  where  A  is  an  amplitude 
factor.  The  correlation  operation  now  becomes 

H(T)«Jg(t)  c(t+T)dt  =/^]_(t)  +  Ag2(t-tQ)J  c(t+T)dt 
=  Jgl(t)  c  (t+T>dt  +  A  Jg2(t-t0)c(t+T)dt-H1(t)  + 

AH2(t-tQ). 

If  the  dispersion  characteristics  of  g-^(t)  and  g2(t)  are 
similar,  as  is  often  the  case  when  their  travel  paths 
are  nearly  the  same,  the  second  peak  will  be  a  replica 
of  the  first  only  delayed  in  time  by  t^  and  modified  in 
amplitude  by  A. 

Theoretically,  H(T)  endures  from  to  +  eo  . 

In  practice,  for  band  limited  data  (15-70  sec  periods, 
for  example),  the  energy  is  nearly  all  contained  in 
t  100  sec  from  the  peak.  A  200  sec  separation  in  the 
arrivals  times  of  different  signals  would  therefore 
result  in  noninterfering  peaks  and  would  allow  the  true 
dispersion  curve  to  be  obtained  by  maximizing  the  real 
part  of  the  transform  of  the  first  peak.  As  time 
separation  decreases  between  the  primary  and  the 


multipathed  arrivals,  the  peaks  will  begin  to  interfere  with 
each  other  and  resolution  will  be  lost.  Empirically  the 
minimum  time  separation  for  which  this  method  yields  accurate 
results  will  be  shown  to  be  less  than  50  sec  for  data  con' 
taining  periods  up  to  70. 


Examples  Using  Synthetic  Signals 

The  relations  just  discussed  can  be  illustrated  empir¬ 
ically  by  the  use  of  synthetic  signals.  Since  paths  from 
Central  Asia  to  the  United  States  are  of  special  interest, 
a  signal  with  dispersion  characteristic  of  such  a  path  was 
chosen  (propagation  distance  approximately=d2, 000  kms).  This 
signal  (Figure  1)  has  a  constant  amplitude  and  a  smoothly 
varying  dispersion  curve  shown  in  Figure  2.  This  dispersion 
curve  (referred  to  as  #1)  was  fixed  at  its  long  period  end 
(3.925  kms/sec  at  75.5  sec)  and  rotated  about  that  point. 

New  curves  were  manufactured  with  velocity  increments  of  .01 
km/sec  at  a  period  of  17  sec.  Six  additional  synthetic 
curves  formed  in  this  way  are  labeled  #2,4,6  on  che  low  ■velo¬ 
city  side  of  #1,  and  3,5,7  on  the  high  side.  A  chirp  wave  form 
was  derived  from  each  of  these  seven  curves  using  the  method 


of  McDonald  et  al  (1974)  .  Each  chirp  was  then  cross-correlated 
with  the  synthetic  signal  to  obtain  seven  correlation  functions 
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with  correponding  numbers.  Correlation  functions  1,  2, 
and  3  are  shown  in  Figures  3-5.  Since  #1  is  a  true 
autocorrelation  function  it  displays  the  characteristic 
of  e.veness  about  the  peak  and  has  the  highest  peak 
magnitude  for  any  of  the  seven  correlations.  See  the 
first  half  of  Table  1,  where  the  magnitude  of  the  peak 
for  #1  was  arbitrarily  assigned  the  value  of  1. 

It  may  also  be  noted  that  correlations  2  and  3  are 
mirror  images  of  each  other  about  the  peak.  This 
characteristic  results  from  the  development  of  a  small 
odd  portion  in  each  of  these  time  series .  Because  chirps 
2-7  have  phase  spectra  unequal  to  the  signal  phase,  the 
imaginary  part  of  the  transform  does  not  vanish  for  these 
correlations  as  it  does  for  #1. 

The  real  and  imaginary  parts  of  the  transform  of 
each  of  the  correlation  functions  were  computed  placing 
tne  maximum  value  of  the  main  peak  at  zero  time  and 
centering  a  200  sec  window  on  this  peak,  then  using  the 
method  of  Filon  (1929)  for  discrete  transforms.  The  real 
part  was  integrated  over  the  known  bandwidth  (72-18  sec 
period-)  for  each  correlation  function.  These  results 
are  listed  in  Table  1  along  with  the  peak  magnitudes, 
and  are  shown  graphically  in  Figure  6.  For  this  particular 
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signal  it  is  clear  that  the  true  dispersion  curve  could 
be  determined  from  either  correlation  function  peak 
magnitude  or  from  the  integrated  real  part  of  the  trans¬ 
formed  correlation  functions.  These  criteria  are  nearly 
equivalent  in  sensitivity. 

An  application  of  the  technique  was  then  illustrated 
using  a  synthetic  signal  with  a  simulated  multipath 
(Figure  7).  This  signal  was  created  nv  delaying  the 
first  synthetic  signal  by  50  sec,  dividing  its  amplitude 
by  2,  zeroing  the  first  100  samples,  and  adding  the 
delayed  portion  to  the  original  signal.  The  true  dis¬ 
persion  curve  of  this  multipathed  signal  is  then  the  same 
as  that  for  the  undisturbed  signal.  The  same  seven 
filters  were  then  applied  to  the  new  signal.  Correlations 
1,  2,  and  3  (Figures  8-10)  show  both  a  main  peak  and  a 
smaller  one  delayed  50  sec,  the  exact  time  shift  of  the 
multipathed  energy.  The  eveness  of  #1  has  been  only 
slightly  disturbed  out  to  the  point  of  the  second  peak. 

An  examination  of  the  sensitivity  (Table  1,  second  part, 
and  Figure  6)  for  this  signal  shows  nearly  the  same 
results  as  were  obtained  from,  the  undisturbed  signal. 

Thus  the  theory  has  been  verified  in  these  synthetic 
examples;  it  has  been  shown  that  the  true  dispersion  curve 
is  indeed  attainable  in  the  presence  of  multipathing . 
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APPLICATION  TO  SEISMIC  DATA 

The  Large  Aperture  Seismic  Array  (LASA)  is  loca<--d 
near  Billings,  Montana.  It  consists  of  21  subarrays  for 
detection  of  high  frequency  seismic  energy  (pattern  shown 
in  Figure  11).  All  of  these  stations  except  those  of  the 
B  ring  include  long  period  seismometers  which  record 
Rayleigh  and  Love  waves.  The  response  curve  for  long 
period  data  is  shown  in  Figure  12.  It  should  be  noted 
that  all  LASA  seismograns  shown  here  are  the  raw  data.  No 
correction  has  been  made  for  instrument  response  nor  has 
there  been  any  prefiltering  of  these  signals. 

The  size  of  this  array  can  cause  serious  problems 
in  the  estimation  of  dispersion  curves  because  multi- 
pathed  signals  often  create  different  effects  at  sensors 
separated  by  more  than  a  few  tens  of  kilometers.  As  a 
result  it  may  appear  that  each  site  has  its  own  particular 
dispersion  curve  even  when  this  is  not  the  case.  For 
example,  E2.  is  only  70  kms  east  of  AO,  thus  if  Rayleigh 
waves  from  an  Asian  earthquake  with  epicenter  10,000  kms 
directly  north  of  the  array  are  recorded  at  both  sites, 
similar  waveforms  would  be  expected.  Figures  13  and  14 


show  the  vertical  channel  data  recorded  at  AO  and  E2 
from  such  an  earthquake  (epicenter  in  Tadzhik,  USSR) . 

An  examination  of  these  two  signals  reveals  little  differ¬ 
ence-  in  the  long  period  Rayleigh  waves;  however,  at 
periods  of  approximately  25  sec  and  less  the  seismograms 
are  quite  different.  A  comparison  of  the  recordings  a+: 

FI  (Figure  15)  and  F3  (Figure  16)  from  this  same  earth¬ 
quake  also  show  variations  from  the  AO  signal.  Either 
structural  differences  across  LASA  or  the  very  small 
differences  in  travel  path  have  caused  various  multipathed 
signals  to  interfere  differently  with  the  first  arriving, 
fundamental  mode  Rayleigh  wave.  Both  the  amplitude  and 
phase  have  been  distorted.  Further  evidence  that  multi¬ 
pathed  energy  is  present  can  be  seen  from  Figures  17  and 
18,  where  the  same  event  is  shown  as  recorded  at  D3 . 

Figure  17  is  the  radial  component  of  motion  at  D3  and 
Figure  18  is  the  transverse.  Since  a  Rayleigh  wave  has 
no  transverse  motion  and  Love  waves  would  have  arrived 
much  sooner  due  to  their  faster  velocity,  the  signal  which 
appears  in  Figure  18  must  have  traveled  other  than  a  great 
circle  path. 

For  the  Tadzhik  earthquake,  the  signals  recorded 
at  AO,  C3,  C4,  Dl,  D2,  D3,  and  D4  are  very  similar  in 


waveform  and  could  be  beam-summed  to  give  a  good  increase 
in  signal  to  noise  ratio  ,  However,  recordings  at  the  E 
and  F  sites  which  constitute  the  outer  rings  at  LASA  show 
significant  effects  of  the  multipathing  problem  which 
would  seriously  degrade  the  beam-sum. 

The  first  step  in  determining  a  dispersion  curve 
for  the  path  from  Tadzhik  to  LASA  was  to  obtain  an  initial 
estimate  of  the  curve.  Figure  19  is  a  display  of  the  AO 
group  velocity  data  computed  from  zero  crossings.  The 
scatter  is  typical  of  real  signals  and  increases  the 
difficulty  in  choosing  a  trial  curve.  A  trial  line  was 
drawn  through  the  data  as  a  starting  point  (broken  line  in 
Figure  19). 

Dispersion  is  clearly  present  in  the  period  range 
of  42  -  16  sec  and  visibly  extends  to  50  and  12  for  some 
sensors  in  the  array.  It  was  determined  that  all  trial 
curves  should  ranga  at  least  to  70  sec  so  that  no  hidden 
long  period  arrivals  would  be  missed.  The  lower  limit 
was  set  at  12  sec.  A  chirp  filter  was  generated  from 
the  first  trial  curve  and  applied  to  the  AO  signal. 

The  resulting  correlation  function  was  then  transformed 
using  the  Filon  method  with  a  200  sec  window  centered  on 
the  maximum  positive  peak,  and  the  real  spectrum  was 
examined  and  integrated  over  the  period  range  of  4.2-16  sec. 


shorter  period  portions  until  an  accurate  estimate  was 
found  for  the  entire  curve.  The  belief  in  the  correctness 
of  this  result  is  reinforced  by  the  even  property  displayed 
in  this  correlation  function  about  its  peak  (see  Figure  20) . 
The  final  curve  is  labeled  7B  and  is  plotted  beside  the 
first  estimate  in  Figure  19. 
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A  sensitivity  diagram  was  constructed  similar  to 
those  for  the  synthetic  sic,;  als  shown  earlier.  With  the 
long  period  end  fixed  the  rest  of  the  curve  was  rotated 
to  obtain  increments  of  velocity  deviation  at  the  short 
period  end.  A  set  of  seven  new  curves  was  created  for 
comparison  with  7B.  The  results  of  integrated  real  spectra 
of  the  correlation  function  transforms  are  listed  in  Table 
2  along  with  correlation  peak  magnitudes  obtained  from 
these  curves.  For  curves  of  faster  velocity  the  sensitivity 
for  this  signal  is  similar  to  that  found  in  the  synthetic 
signals,  i.e.  relative  magnitudes  drop  as  the  velocity 
deviation  is  increased  at  the  short  period  end.  For  the 
slower  velocity  curves,  however,  a  secondary  maximum 
exists  below  the  correct  curve.  At  a  short  period  devia¬ 
tion  of  -.03  km/sec  the  relative  magnitude  rises  to  .93 
after  an  initial  drop  to  .85.  This  phenomenom  can  most 
likely  be  explained  by  noting  that  as  the  curve  is 
rotated  from  its  correct  value  in  Figure  19  it  first  goes 
through  a  space  where  little  dispersion  data  exists;  as 
rotation  is  continued,  a  cluster  of  low  velocity  points 
is  encountered  which  increases  the  correlation  function 
peak  magnitude.  Having  passed  this  area  the  magnitudes 
fall  once  again.  See  Figure  21. 


Using  7B  as  a  starting  point,  the  dispersion 
curve  was  then  determined  for  other  sensors  in  the  LASA 
array.  It  was  found  that  ?n  ^as  also  the  most  accurate 
estimate  for  all  of  the  C  and  D  sites  as  well  as  El,  E3, 
and  FI.  A  slightly  different  dispersion  curve  (labeled 
17B )  was  found  for  E2,  F2,  and  F3  and  a  third  variation 
of  7B  (18B)  was  judged  the  best  estimate  for  E4  and  F4. 

The  three  curves  7B,  17B,  and  18B  are  very  similar. 
The  only  differences  occur  in  the  period  range  29-34  sec. 
These  variations  are  magnified  for  comparison  in  B'igure 
22.  A  complete  list  of  the  integrated  real  parts  of  the 
transformed  correlation  functions  are  given  in  Table  3  for 
these  three  curves  applied  over  the  entire  array.  A  geo¬ 
graphical  relationship  between  the  sites  and  their  corre¬ 
sponding  dispersion  curves  is  shown  in  Figure  23. 

One  of  the  secondary  goals  of  this  study  was  to 
make  use  of  the  dispersion  curves  determined  from  this 
technique  by  creating  a  chirp  filter  which  would  be  useful 
in  amplifying  signal  to  noise  ratio  for  small  magnitude 
events  near  Tadzhik.  It  is  suggested  that  any  of  these 
curves —  7B,  17B,  or  18B  — should  give  a  useful  filter 


for  this  purpose. 
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This  idea  was  tested  using  a  small  magnitude  earth¬ 
quake  whose  epicenter  was  located  approximately  1400  kms 
northeast  of  the  Tadzhik  event.  The  signal  in  Figure  24 
is  barely  recognizable  on  the  AO  vertical  channel.  Figure 
25  is  the  resulting  correlation  function  with  filter  18B 
applied  to  the  signal.  The  increase  in  signal  to  noise 
ratio  is  apparent.  Both  the  least  time  arrival  and  a 
multipath  or  after  shock  delayed  140  sec  are  now  clearly 
visible . 

The  one  difficulty  encountered  by  this  method  in 
its  present  form  is  the  amount  of  labor  required  in  com¬ 
paring  all  the  trial  curves  brought  about  by  the  tremen¬ 
dous  number  of  perturbations  possible  for  any  initial 
estimate.  An  automated  system  of  curve  variation  from  a 
preliminary  estimate  is  needed.  One  way  of  satisfying 
this  requirement  is  therefore  suggested.  If  the  initial 
curve  is  modified  to  a  time  function  dependent  on 
frequency  with  a  fixed  intercept  on  the  time  axis,  it 
would  then  be  possible  to  compute  the  coefficients  of  a 
fifth  or  sixth  order  polynomial  which  would  closely 
approximate  the  curve.  Let  U(f)  =  a+bf+cf ^+df 3+ . . . 

'-htn  the  velocity  intercept  (=a)  could  be  fixed  while  the 


initial  slope  (b)  and  the  coefficient  which  determines 
the  inflection  point  (c)  would  be  varied  slightly — 
enough  to  perturb  the  high  frequency  end  of  the  curve 
through  a  range  of  ±  .03  or  .04  kms/sec.  It  is  possible 
that  a  set  of  curves  generated  in  this  manner  could  have 
their  chirps  and  signal  correlations  examined  quickly  and 
efficiently  by  a  computer  which  would  then  search  for  the 
largest  peak  magnitude  among  the  correlations  and  thus 
determine  the  best  curve  in  one  complete  operation. 
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CQNCLUS IONS 

It  has  been  clearly  demonstrated  that  multipathing 
causes  distortions  in  both  the  amplitude  and  phase  of 
dispersed  wavetrains.  The  technique  proposed  here  for 
determining  dispersion  curves  has  been  shown  through  the 
use  of  synthetic  signals  to  give  nearly  identical  results 
for  both  undisturbed  and  multipathed  data.  In  both  cases 
the  trial  curves  when  fixed  at  the  long  period  end,  could 
be  used  to  determine  the  correct  curve  to  within  .01  km/sec 
deviation  at  the  short  period  end.  Also,  evidence  indicates 
t>_,t  these  results  were  duplicated  when  the  method  was 
applied  to  LASA  Rayleigh  waves  in  the  actual  determination 
of  an  unknown  dispersion  curve. 

As  demonstrated,  the  comparative  accuracy  of  this 
technique  far  exceeds  that  expected  from  other  methods 
when  multipathed  interference  is  present. 
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Relative  magnitudes  of  correlation  functions  for  undisturbed 
synthetic  signal 
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Relative  magnitudes  of  correlation  functions  for  multipathed 
synthetic  signal 
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Table  2 


Relative  magnitudes  of  correlation  functions  for  AO 
Tadzhik  signal 


Short  Period  Dev. 

from  Curve  7B  Max.  Peak  Real  Integral 
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+  .02 
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TABLE  3 
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WEAL  INTEGRAL  0F  C9WRELATI0N  FUNCTION  FOR  TADZHIK  EARTHq JAj 
OF  17  MARCH  1972  9-17-10*  .  1 


SITE 
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D3 

176 

1»17 

•j 

D3 

18U 
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04 
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1  •  44 

1 
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24 
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2 
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4 
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Synthetic  signal.  1 

■> 

Dispersion  curve  (#1)  of  synthetic  signal  1 

shown  in  Figure  1 .  | 

Correlation  function  of  synthetic  signal  I 

with  chirp  1.  (true  autocorrelation)  j 

Correlation  function  of  synthetic  signal  k 

with  chirp  2 .  | 

Correlation  function  of  synthetic  signal  with  | 

chirp  3 . 

4 

I 

Sensitivity  plot.  Relative  magnitudes  of  j 

the  correlation  function  peaks  and  the  real  j 

spectral  integrals  for  the  seven  chirps  ] 

applied  to  (1)  the  undisturbed  synthetic  \ 

signal  and  (2)  the  multipathed  synthetic  \ 

signal.  I 
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9 

10 

11 

12 


Multipathed  signal.  The  synthetic  signal  had 
its  first  100  seconds  zeroed;  it  was  then 
shifted  50  seconds  and  added  at  half  ampli¬ 
tude  to  the  original  signal. 

Correlation  function  of  multipathed  signal 
with  chirp  1. 

Correlation  function  of  multipathed  signal 
with  chirp  2. 

Correlation  function  of  multipathed  signal 
with  chirp  3 . 

LASA  array  map. 

LASA  long  period  response  curve. 
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13  Tadzhik  earthquake — AO  vertical  channel.  4 


14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 


Tadzhik  earthquake — E2  vertical  channel. 

Tadzhik  earthquake — FI  vertical  channel. 

Tadzhik  earthquake — F3  vertical  c1  .inel. 

Tadzhik  earthquake — D3  radial  channel. 

Tadzhik  earthquake — D3  transverse  channel. 

AO  Tadzhik  group  velocity  plot.  The  data 
points  are  computed  from  zero  crossings  of 
the  signal.  Curve  7B  shows  its  relation 
to  the  dispersion  data. 

.Cross  correlation  function  of  AO  Tadzhik 
vertical  signal  and  chirp  7B . 

Sensitivity  plot.  Relative  magnitudes  of 
correlation  function  peaks. 

Curves  7B,  17B,  and  18B .  The  differences 
between  these  curves  has  been  accentuated 
by  the  large  scale  plot. 

Sensitivity  plot.  Relative  magnitudes  of 
the  correlation  function  peaks  and  the  real 
spectral  integrals  for  the  eight  chirps 
applied  to  the  AO  Tadzhik  vertical  signal. 

AO  vertical  channel  for  an  earthquake  approx¬ 
imately  1400  kms  northeast  of  the  Tadzhik 
event . 

Correlation  function  of  event  from  Figure  25 
with  filter  18B  applied. 
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